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1. Background and Motivation 



Abstract. We construct a Lax pair for the g-Painleve system from first principles by employing the 
general theory of semi-classical orthogonal polynomial systems characterised by divided-difference operators 
on discrete, quadratic lattices |16|. Our study treats one special case of such lattices - the (j-lrnear lattice - 
through a natural generalisation of the big g-Jacobi weight. As a by-product of our construction we derive 
the coupled first order q-difference equations for the Eg <j-Painleve system, thus verifying our identification. 
Finally we establish the correspondences of our result with the Lax pairs given earlier and separately by Sakai 
. and Yamada, through explicit transformations. 

o 

C; 

Since the recent discoveries of ^-analogues of the Painleve equations, see for example [H and Ifl2l 
which are of relevance to the present study, and their classification (of these and others) according to 
the theory of rational surfaces by Sakai [13 1 interest has grown in finding Lax pairs for these systems. 
This problem also has the independent interest as a search for discrete and ^-analogues to the isomon- 
odromic systems of the continuous Painleve equations, and an appropriate analogue to the concept of 
monodromy Such interest, in fact, goes back to the period when the discrete analogues of the Painleve 
. equations were first discussed, as one can see in fTTl . 

In this work we illustrate a general method for constructing Lax pairs for all the systems in the Sakai 
scheme, as given in the study 1161 , with the particular case of the E 6 system. In this method all aspects 
of the Lax pairs are constructed, and in the end we verify the identification with the Eg system by 

deriving the appropriate coupled first order ^-difference equations. We will utilise the form of the Eg 
^-Painleve system as given in 1 5 1 and [ 6 1 in terms of the variables /, g under the mapping 

(t,f,g)»(qt, f(qt)=f,g(qt) = g), 



and f(q l t) = f, etc. In these variables the coupled first order ^-difference equations are 

(i i) (gf- i)(gf - 1) = ^s-Wis-Dihg-Wig-i) 



(Mi (1-2) (fi-l)(fg-l) = qi 



,r* ^(/-W- b 2 )(f-h)(f-h) 



(f-b 5 qt)(f-b-H) 

with five independent parameters b\, . . . , b(, subject to the constraint b^b^bi - 1. 

Our approach is to construct explicitly a sequence of T-functions starting with a deformation of a 
specific weight in the Askey table of hypergeometric orthogonal polynomial systems Q. The weight 
that we will take is the big ^-Jacobi weight given by Eq. (14.5.2) of (7| 

(1.3) w(x) = — . 

(x,bc i x;q) co 

The essential property of this weight, and the others in the Askey table, that we will utilise is that they 
possess the ^-analogue of the semi-classical property with respect to x, namely that it satisfies the linear, 
first-order homogeneous (j-difference equation 

wjqx) _ a(l - x)(c - bx) 
zv(x) (a - x)(c - x) 

where the right-hand side is manifestly rational in x. Another feature of this weight is that the discrete 
lattice forming the support for the orthogonal polynomial system is the (j-linear lattice, one of four 
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discrete quadratic lattices. Consequently the perspective provided by our theoretical approach, then 
indicates that this case is the master case for the ^-linear lattices (as opposed to the d! system, for 



example) and all systems with such support will be degenerations of it. The weight il.3i has to be 
generalised, or deformed, in order to become relevant to ^-Painleve systems, and such a generalisation 
turns out to introduce a new variable t and associated parameter so that it retains the semi-classical 
character with respect to this variable. Using such a sequence of T-f unctions one employs arguments to 
construct three systems of linear divided-difference equations which in turn characterise these. One of 
these is the three-term recurrence relation of the polynomials orthogonal with respect to the deformed 
weight, which in the Painleve theory context is a distinguished Schlesinger transformation, while the 
two others are our Lax pairs with respect to the spectral variable x and the deformation variable t. Our 
method constructs a specific sequence of classical solutions to the E 6 system and thus is technically 
valid for integer values of a particular parameter, however we can simply analytically continue our 
results to the general case. 

Lax pairs have been found for the E?' system system using completely different techniques. In IIT51 
Sakai used a particular degeneration of a two-variable Gamier extension to the Lax pairs for the D^ 1 ' 
(^-Painleve systerrQ (see [14 [ for details on the multi-variable Gamier extension). More recently Yamada 
|[T7l has reported Lax pairs for the E^ system by employing a degeneration starting from a Lax pair for 

the Eg 1 ' ^-Painleve equation through a sequence of limits E g ^ — » Ey* — > E^ 1 '. 

The plan of our study is as follows. In Section|2]we recount the notations, definitions and basic facts 
of the general theory Ifl6l in a self-contained manner omitting proofs. We draw heavily upon this theory 
in Section [3] where we apply it to the ^-linear lattice and a natural extension or deformation of the big 
q-Jacobi weight. Again, using techniques first expounded in [16 1, we find explicit forms for the Lax 
pairs and verify the identification with the Ej. ^-Painleve system. At the conclusion of our study, in 
SectionHJ we relate our Lax pairs with those of both Sakai and Yamada. 

2. Deformed Semi-classical OPS on Quadratic Lattices 

We begin by summarising the key results of [16J, in particular Sections 2, 3, 4 and 6 of that work, which 
relate to semi-classical orthogonal polynomial systems with support on discrete, quadratic lattices. In 
the divided difference calculus of discrete non-uniform lattices of quadratic type [9[ the divided-difference 
operator (DDO) D x is defined by 

(Z1) = L+ (x)-u(x) ' 

where the advanced and retarded nodes i±(x) are characterised by the condition that D Y : U„[x] — » 
n„_i[x], V n > 0. Here n„[x] is the linear space of polynomials with degree at most n over x e C. This 
implies that i±(x) are the y-roots of of the quadratic equation 

(2.2) J{y 2 + 2Sxy + Cx 2 + 2<Dy + 26x + T = 0. 

For the quadratic lattices we are also obliged to define an average operator M v by 

(2.3) M % /w = i[/a + w)+/a-W)]- 

The dual nodes and their difference are also denoted by 

(2.4) y ± = i ± (x), 

(2.5) Ay(x) = i+(x) - u(x). 

The quadratic lattices are classified into four sub-cases (9j: ^-quadratic, quadratic, ^-linear and linear, 
as the conic sections are divided into the elliptic/hyperbolic, parabolic, intersecting straight lines and 
parallel straight lines respectively. We define the ID-Integral of a function defined on a quadratic x-lattice 
G by f : G — » C, whether finite or denumerable, by the Riemann sum over the lattice points 

(2.6) f Dx/(x) = y Ay(x)f(x). 

All of the results we give in this section will apply to any of the quadratic lattices enumerated above. 



this later system is also known as the <J-Pvi system and its Lax pairs were constructed in 1 4 [ 
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Central to our study are orthogonal polynomial systems (OPS) defined on G, and a general reference 
for a background on these and other considerations is the monograph by Ismail |3[. Our OPS is defined 
via orthogonality relations with support on G 



(2.7) 



J" Dx w(x)p„(x) l m {x) 



J 0, < m < n 
\h n , m = n 



n>0, hn*0, 



where [lm{x)}^_ is arw system of polynomial bases with exact deg v (/,„) = m. Such relations define 
a sequence of orthogonal polynomials {p>i(*)}^l under suitable conditions (see [3|). An immediate 
consequence of orthogonality is that the orthogonal polynomials satisfy a three term recurrence relation 
of the form 



(2.8) 



a n +ip n +i(x) = (x- b n )p n (x) - a n p n -i(x), n>0, a„£0, p-\ = 0, p = y . 



However we require non-polynomial solutions to this linear second-order difference equation, which 
are linearly independent of the polynomial solutions. To this end we define the Stieltjes function 



C w(y) 
f(x)= Dy^^, x?G. 



(2.9) 

A set of non-polynomial solutions to j2.8b , termed associated functions or functions of the second kind, and 
which generalise the Stieltjes function, are given by 



(2.10) 



ln(x) = 



Dy w(y) : 



y 



n>0, x<£G. 



With the polynomial and non-polynomial solutions we form the 2x2 matrix variable, which occupies 
a primary position in our theory: 

q n (x) 

Fnw 

(2.11) 



Y n (x) 



7„-l(x) 



w{x) 

q n -i(x) 



w(x) 



In this matrix variable the three-term recurrence relation takes the form 

1 (x — b n —a„ 
a n+ i \ On+l 



(2.12) 



K n (x) 







de tX„ = -US-. 



A key result required in the analysis of OPS are the expansions of polynomial solutions about the fixed 
singularity at x - oo 



(2.13) 



Pn(x) = y„ 



x" - 



(n-\ \ 
\i=0 



x"- 1 + 



11-1 



0<i<j<n 



v n-2 



+ 0(x n " 3 ) 



valid for n>l, while for the associated functions the expansions read 



(2.14) 



<7n(*) = Jn 



, i'=0 



x-"- 2 + 



11+1 



^0<t<;'<n 



-n-3 



+ 0(x""- 4 ) 



valid for n > 0. 

In order to proceed any further we need to impose some structure on the weight characterising our 
OPS — in particular its spectral characteristics — and this takes the form of the definition of a D-semz- 
classical weight [9|. Such a weight satisfies a first-order homogeneous divided-difference equation 

w(y+) W + Ay V 



(2.15) 



w(y_) W-AyV 

where W(x), V(x) are irreducible polynomials in the spectral variable x. As a consequence of this, under 
reasonable assumptions on parameters of the weight, the Stieltjes function satisfies an inhomogeneous 
form of H2.15} 

(2.16) WD x f = 2VM x f + U, 

where in addition U(x) is a polynomial of x. A generic or regular D-semi-classical weight has two properties: 

(i) strict inequalities in the degrees of the spectral data polynomials, i.e., deg W = M, deg V — M—l 
and deg^U = M — 2, and 
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(ii) the lattice generated by any zero of ( W 2 - Ay 2 V 2 )(x) does not coincide with the lattice generated 
by another zero, i.e., for any x such that (W 2 - Ay 2 V 2 )(x) = then (W 2 - Ay 2 V 2 )(i 2Z x) + 0. 

Further consequences of semi-classical assumptions are a system of spectral divided-difference equa- 
tions for the matrix variable Y n , i.e., the spectral divided-difference equation 

(2.17) D x Y n (x) := A n M x Y n (x) = -J— ( -"f n{ $ 7l \m x Y n (x), n > 0, 



W n (x) V?„e„_i(x) -Q„(x) - 2V(x) 

with A„ termed the spectral matrix. For the D-semi-classical class of weights the coefficients appearing 
in the spectral matrix, W n , Q„, ©„, are polynomials in x, with fixed degrees independent of the index n. 
The spectral coefficients have terminating expansions about x = oo with the leading order terms 

(2.18) W„(x) = ±W+ i[W+ AyV]|^J + \[W - AyV] + 0{x M ~ l ), n > 0, 

(2.19) @ n {x) = [W + AyV]l^] ^— [W - AyV] (— ) + 0(x M " 3 ), n > 0, 

(2.20) n n{x ) + v(x) = ^[W + AyV]^- ) j - ±-[ W - AyV](^-\ +0(x M " 2 ), n > 0, 
where M = deg v (VV„). 

Compatibility of the spectral divided-difference equations J2.17II and recurrence relations (|2.12b imply 
that the spectral matrix and the recurrence matrix satisfy 

(2.21) K n (y + ) (l - \ Ay (l + \AyA n ) = (l- \Ay A n+l )~ l (l + \AyA n+l ) K n (yJ), n > 0. 

These relations can be rewritten in terms of the spectral coefficients arising in l!2.17b as recurrence 
relations in n, 

(2.22) W„ + i = W„ + \Ay 2 ®„, n>0, 

(2.23) Q„ + i + Q„ + 2V = (M x x - b„)®„, n>0, 
(2.24) 

(WnQn+i - W„ +1 Q„)(M,x - b n ) = -iAy 2 Q )1+1 Q H + W„W n+1 + a 2 n+l W„0 n+1 - a 2 n W n+1 ® n -i, n>0. 

Another important deduction from these relations is that the spectral coefficients satisfy a bilinear 
relation 

(2.25) W„(W„-W) = -lA y 2 det( flH g ; _ i n > 0. 

The matrix product appearing in (|2.2H , and recurring subsequently, has the evaluation 

(2.26) (l-\AyA n )~ X (l + \AyA n ) 

1 (2W n -W + Ay(Q n + V) -Aya n &„ \ 

~W + AyV\ Ay a n @ n - X 2W„ - W - Ay (Q„ + V))' 

This result motivates the following definitions 

(2.27) 2B ± := 2W„ - W ± Ay (Q„ + V), n>l, 

(2.28) %+ := Aya„®„, %--.=■ Ay a n © n -\, n>l, 

whilst for n — we have 2B ± (n = 0) := W ± Ay V, % + {n = 0) := -Aya ylU, and X-(n = 0) := 0. Thus we 
define 



(2.29) A; := 



9B + -X, 



In a scalar formulation of the matrix linear divided-difference equation (|2.17b one of the components, p n 
say, satisfies a linear second-order divided-difference equation of the form 

W + AyV\ , (W-AyV\ , 



(2.30) 



11 Aye,, 



(t+(x)) + 



2B_ 



Aye,, 



(u(x)) \p n (i + o t _(x)) = 0. 



Thus far our theoretical construction can only account for the OPS occurring in the Askey table 
- the hypergeometric and basic hypergeometric orthogonal polynomial systems Q. To step beyond 
these, and in particular to make contact with the discrete Painleve systems, one has to introduce pairs 
of deformation variables and parameters into the OPS. We denote such a single deformation variable 
by t, defined on a quadratic lattice (and possibly distinct from that of the spectral variable), with 
advanced and retarded nodes at t±(f) = u+, Au = i+(t) — i-(t). We introduce such deformations with 
imposed structures that are analogous to those of the spectral variable. Thus, corresponding to the 
definition l!2.15l l, we deem that a deformed D-semi-classical weight w(x; f) satisfies the additional first-order 
homogeneous divided-difference equation 

(2.3D ^ = «±£ffc t) , 



w{pc; U- 



R-AuS 



where the deformation data polynomials, R(x; t), S(x; t), are irreducible polynomials in x. The spectral 
data polynomials, W(x; t), V(x; t), and the deformation data polynomials, R(x; t), S(x; f), now must satisfy 
the compatibility relation 



(2.32) 



W + AyV, R + AuS, , W + AyV R + AuS, 

;(*;«+) D — (y-;*) = w —, (x;u-) j—^ (y + -,t). 



W-AyV" ' ^R-AuS" J ' W - AyV 

The deformed D-semi-classical deformation condition that corresponds to l|2.161 is that the Stieltjes 
transform satisfies an inhomogeneous version of i2.31|l 



(2.33) 



RD t f = 2SM t f + T, 



with T(x; f) being an irreducible polynomial in x with respect to R and S. Compatibility of spectral and 
deformation divided-difference equations for / implies the following identity on IT and T 



(2.34) Ay 



(W + AyV)(x; u+) 



(W + AyV){x;u-) 



(R + A«S)(y_; t)U(x; u_) - (R - A«S)(y_; t)U(x; u+) 



Au 



( W + AyV)(x- U+ )T(y_; f) - (W - AyV){x; «+) 



(R- AuS){y-;t) 



T{y + ;t) 



(R - AuS)(y + ; t) 

Corresponding to the Il2.17b the deformed D-semi-classical OPS satisfies the deformation divided-difference 
equation 

<$>n\ 



(2.35) 



D f Y„ := B n M t Y n = 1. IE? _ 7 



M t Y„, n > 0. 



The deformation coefficients appearing in matrix B„ above satisfy a linear identity 

a„ 



(2.36) 
and a trace identity 

(2.37) 



W« = 



Au (P, + S„) = 2H„ 



n > 1, 



a n -\ 



R + AuS R - AuS 



a n (u-) a„(u+) 



n > 0, 



which means that only three of these are independent. Here H n is a constant with respect to x and 
arises as a decoupling constant which will be set subsequently in applications to a convenient value. 
The deformation coefficients are all polynomials in x, with fixed degrees independent of the index n but 
with non-trivial f dependence. Let L = max(deg a .i?, deg^S). As x — > oo we have the leading orders of the 
terminating expansions of the following deformation coefficients 



(2.38) 



TrRn = -(Yn(u+) + Yn{U-)) 



R 



AuS R + AuS 
+ 



y n -i(u+) 
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yn-i(w-) 



+ 0(x L_1 ) n>0, 



(2.39) 
and 
(2.40) 



Ait 

2H~n 



7n(U-) 



7n(u+) 



x' 1 + 0(x L_2 ), n > 0, 



^r„ = ()/„(«_) - y n (u+)) 



R + AuS R - AuS 
+ 



+ 0(x L_1 ), n > 0. 



y«-i(«-) y«-i(w+) 

Compatibility of the deformation divided-difference equation l|2.35b and the recurrence relation l|2.12t 
implies the relation 

(2.41) K n (; u + ) (l - \ AuBn)' 1 (l + |AmB„) = (l - \AuB n+1 )~ l (l + I AuB n+1 ) K„(; «_), n>0. 

From this we can deduce that the deformation coefficients, R n , T n , O n , satisfy recurrence relations in n in 
parallel to those of (|2.22l2.23b 

cirn--\(u-) a„(u-) 

(2.42) "+ n \ -2R n+l + Aur„ +1 ) + -^-^(2R„ + AuT n ) 



H„ 



+i 



r , / ^ A »^ n / J R + AuS R - AuS 

-[x -b n (m_ ]— O n + 2fl n (M_) — — — - 

H„ \ fl„(M_) «„(«+) 



n > 0, 



(2.43) ^ll(2R H+1 + Aur n+1 ) + ^(-2R„ + AuT„) 



r , / m Am ^ n / jR + AuS R-AuS 

= -[x- b„(u + )]—<J? n + 2a n (u + ) — — — - 

H n \ a n (u-) a n (u+) 

The deformation coefficients satisfy the bilinear or determinantal identity 



n > 0. 



(2.44) 



Rl + \Au 2 [r„H„ - ® n W n ] = -H„R„ 



R + AwS R - AuS 
+ 



, «>0, 



a„(w_) a„(u+) 

which is the analogue of J2.25II . The matrix product given in j2.41|l has the evaluation 
(2.45) (l - \Au B,,)" 1 (l + \Au B n ) 

a„(u-) 



( R- AuS 
2R„ + 2H n : — — + Au T„ 



2H n (R + AuS) 

This again motivates the definitions 
(2.46) 



a n {u+) 
AuW n 



Au O n 

R + AuS . „ 
2R„ + 2H„ — — - Am r„. 



a n {u-) 



n>0. 



(2.47) 
(2.48) 

valid for « > 1 together with 
(2.49) 



5R+ := 2R„ + 2H n R + Au T nr 

a n {u+) 

91_ := 2R„ + 2H n R + AuS - Au Y n , 
a„(u-) 

^ + :=-AuO„, *P_:=AuW„, 



b: := 



9U -^P + 



Our final relation expresses the compatibility of the spectral and deformation divided-difference 
equations. The spectral matrix A n (x;t) and the deformation matrix B n (x;t) satisfy the D-Schlesinger 
equation 

(2.50) (l - \ AxjA n {; M+ )) _1 (l + \AyA n {; u + )) (l - \AuB n {y.; ))"' (l + \AuB n {y.; )) 

= (l - jAuB n (y + ; ))"' (l + ±A«B„(y + ; )) (l - \AyA n {; w_)) _1 (l + \AyA n (; «_)) . 
Let us define the quotient 

(W + Ay V)(x; «+) (R + A«S)(y_; f) (W - AyV)(x; u+) (R - AuS)(y_; f) 



(2.51) 



A" = 



(W + AyV)(x; u_) (R + AuS)(y+; t) (W - AyV)(x; «_) (R - AwS)(y + ; f) ' 



The compatibility relation (12.50ft can be rewritten as the matrix equation 

(2.52) x K(y + ; t)A* n {x; m_) = A* n (x; u + )B* n (y-, t), 
or component-wise with the new variables in the more practical form as 

(2.53) X [28 + (x; U _)3? + (y + ; t) - 2L(z; U-W + (y + ; t)] = 2B + (x; w + )^ + (y-; - Z + {x; M + )^_(y_; t), 

(2.54) X [X+te u-)5R + (y + ; + 2B_(x; U-W + (y + ; t)] = Z + (x; M+ )5R-(y_; f) + 2B + (x; u + )<J3 + (y_; t), 

(2.55) x [3L(*; H-)W_(y + ; f) + m + {x; «_)<P_(y + ; f)] = 2L(*; « + )5R + (y-; + H + )^-(y_; f), 

(2.56) ^ [2IL(x; H_)M_(y + ; t) - 2 + (*; «-MP-(y+; 0] = i*+)SR_(y-; f) - 3L(ac; M + )^ + (y_; f). 

For a general quadratic lattice there exists two fixed points defined by L+(x) = and let us denote 
these two points of the x-lattice by Xi and Xr. By analogy with the linear lattices we conjecture the 
existence of fundamental solutions to the spectral divided-difference equation about x = Xi, Xr which 
we denote by Yi, Yr respectively. Furthermore let us define the connection matrix 

(2.57) P{x;t):=Y R {x;ty l Y L {x;t). 



From the spectral divided-difference equation l|2.17[l it is clear that P is a D-constant function with 
respect to x, that is to say 

(2.58) P(y + ;t) = P(y-;t). 

In addition it is clear from the deformation divided-difference equation l|2.35b that this type of deforma- 
tion is also a connection preserving deformation in the sense that 

(2.59) P(x;u+) = P(x;u-). 

This is our analogue of the monodromy matrix and generalises the connection matrix of Birkhoff and 
his school m, (21/ although we emphasise that we have made an empirical observation of this fact and 
not provided any rigorous statement of it. 

3. Big ^-Jacobi OPS 

As our central reference on the Askey table of basic hypergeometric orthogonal polynomial systems 
we employ |8], or its modern version |7J. We consider a sub-case of the quadratic lattices, in particular 
the ^-linear lattice in both the spectral and deformation variables x and t in its standardised form, so 
that i+(x) = qx, u(x) = x, Ay(x) = (q — l)x and f+(f) = qt, u(t) = t, Au(t) = (q - l)f. In [7] the big rj-Jacobi 
weight given by Eq. (14.5.2) is 

(3.1) w(x) = \ ' - ' H > , 

(x,bc i x;q) OD 

subject to < aq, bq < 1, c < with respect to the Thomae-Jackson ^-integral 

(3.2) dp f{x). 

Jbq 

The (/--shifted factorials have the standard definition 

OO 

(3.3) (fl/'(?)oo = T|(l |^| < 1, («i,...,«n;<?)«i = (di',q)oa-.-ia n ;q) 00 . 

We deform this weight by introducing an extra ^-shifted factorial into the numerator and denominator 
containing the deformation variable and parameter, and relabeling the big ^-Jacobi parameters. We 
propose the following weight 

(b 2 x / b 3 x,b- 1 xr 1 ;q) 0o 
(bix,biX,b 6 xt L ;q) oa 

A condition b^b^b^ = 1 will apply, so we have four free parameters. We do not need to specify the 
support for this weight for the purposes of our work, but suffice it to say that any Thomae-Jackson 
^-integral with terminals coinciding with any pair of zeros and poles of the weight would be suitable. 
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The spectral data polynomials are computed to be 

(3.5) W + Ay V = b 6 (l- b x x) (1 - b 4 x) (t - b 6 x) , 

(3.6) W - AyV = (1 - b 2 x) (1 - b 3 x) (b 6 t - x) . 

Clearly the regular M = 3 case is applicable and we seek solutions to the spectral coefficients with 
deg x W„ = 3, deg v Q„ = 2,deg 3 .0„ - 1. Our procedure is to employ the following algorithm, as detailed 
in 1 16 1. Firstly we parameterise the spectral matrix in a minimal way; secondly we relate the parameter- 
isation of the deformation matrix to that of the spectral matrix and thus close the system of unknowns; 
and finally utilise these parameterisations in the system of over-determined equations to derive evolu- 
tion equations for our primary variables. What constitutes the primary variables will emerge from the 
calculations themselves. 

Proposition 3.1. Let us define a new parameter b$ replacing q n by 

(3.7) f = J^, neZ, . 

Let the parameters satisfy the conditions q + \,b$ + q ll2 ,±l,q 112 , b\b^ + 0, oo and b 2 b 3 + 0, oo. Given the 
degrees of the spectral coefficients we parameterise these by 

(3.8) 2W n - W = w 3 x 3 + w 2 x 2 + w\x + w , 

(3.9) Q„ + V = v 2 x 2 + v\x + vo, 

(3.10) e„ = u x (x - A„) . 

Let A n be the unique zero of the (1,2) component of A* n , i.e., ©„(x) and define the further variables v n = 
(2W„ - W)(A„, t) and \i n = (Q„ + V)(A„, t). Then the spectral matrix is given by 



(3.11) 2W n -W = x 2 ^ + l(x-A„) 



(3.12) G n + V = n n 
b, 



l + (b 1 + b 2 + b 3 + b i )b 6 t + b 2 x + \ n 

% V + : -x - 2tb 6 —-^- 

A n At 



— 6 (x - An) - (1 - b\) 2 A\x - 2b\ [b- 1 + b- 1 + b~ l + bl 1 + (h + bf)t - 2A„] A 2 n 

+ b 5 b-\l + bl) [(l + (h + b 2 + b 3 + h) b 6 t + bl) A n + 2v„ - 2tb 6 ] \, 



and 

b 6 (l - qbl) 

(3.13) ®» = - n ' 5 (x-An). 

q(l - q)b 5 

We note that A„, y. n , v„ satisfy the quadratic relation 

(3.14) v\ = (1 - q) 2 A 2 n u 2 n + b 6 (hA„ - l)(b 2 A„ - l){b 3 A„ - l)(b 4 A n - 1)(A„ - tb 6 )(b 6 A n - t). 

Proof. Consistent with the known data, i.e., the degrees, from (|2. 1812. 1912.20b we compute the leading 
coefficients to be 

h(l - qb 2 ^ 



(3.15) m = - 

(3.16) v 2 = - 

(3.17) w 3 = - 



q(l - q)b 5 ' 
be{l-bl) 
2(1 - # 5 ' 
b 6 (l + bl) 



2b 5 

confirming the relation given by the coefficient of [x 6 ] in (12.25b , — (q — l) 2 v^ + bl . In addition we 
identify the diagonal elements of the [x 3 ] coefficient of A* 

b 

(3.18) k+ = zv 3 + (q - l)v 2 = -b 5 b 6 , k_ = w 3 - (q - l)v 2 = ~. 
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From the coefficient of [x°] in (12.25b we deduce (modulo a sign ambiguity) 

(3.19) w = b 6 t, 
and from the coefficient of [x 1 ] in l|2.25b we similarly find 

(3.20) u>i = -\ [l + t{h + b 2 + b 3 + h)b 6 + bj] . 
Now utilising the condition v„ = (2W„ - W)(A„, t) we invert this to compute 



l + t{b 1 + b 2 + b 3 + b i )b 6 + b 2 v n -tb 6 

(3.21) w 2 = - + \h(h + b^An + b . 

Proceeding further we infer from the coefficient of [x 5 ] in d2.25|) that 

(l + bj) w 2 - (b- 1 + b~ l + b- 1 + b-^hh - b 5 {l + b\)t 

(3.22) Vl = i ^ j \ — , 

(l-?)(l-^) 

and employing the previous result for iv 2 we derive 

(i - b 2 ) b 5 b 6 h- 1 + b- 2 l + b- 1 + bf + (h + b-M 

(3.23) (1-q)) L 1 



l + bl 
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1 + t(h + b 2 + b 3 + b 4 )b 6 + b\ b 6 (l + b\) Vn - tb 6 

+ — ^ A n + 



2A„ 2b 5 A\ ' 

This leaves Vq to be determined. Imposing the relation \i n = (Q„ + V)(A„, t) we can invert this and find 

(3.24) (1 - q) (l - bj) v = (1 - q) (l - bj) u„ + b 5 b 6 [b~ l + fc" 1 + b~ l + b^ 1 + (b h + b?)t] A n - 2b 5 b 6 A 2 n 

- 1(1 + bj) [l + t(h+b 2 + b 3 + b 4 ) h + bj] + ( —-^ (b 6 t - v„) . 

This concludes our proof. □ 

Remark 3.1. We observe that the appearance of the quantity q n b\bnb( ! with n e Z>o and its replacement by 
the new parameter b$ constitutes a special condition. This condition is one of the necessary conditions 
for a member of our particular sequence of classical solutions to the (j-Painleve equations, and 
is built-in by our construction. The other condition derives from the initial conditions n = in our 
construction, see j2.8i and following (|2.28b . 

From our deformed weight 03 .4b we compute the deformation data polynomials to be 

1 

(3.25) R + AuS = — (b 6 qt - x) , 

b(, 

(3.26) R - AuS = (qt - b 6 x) . 

We can verify that the compatibility relation (12.32b is identically satisfied by our spectral and deformation 
data polynomials. We see that this places us in the class L = 1. We will employ an abbreviation for the 
dependent variables evaluated at advanced or retarded times, e.g., 

(3.27) A„(f) = A„, A n (qt)=X n , A n {q- 1 t) = A n . 

In the second stage of our algorithm we parameterise the Cayley transform of the deformation matrix 

(3.28) B n = fr n>0, 

consistent with known degrees, i.e., deg^Si* = 1, deg Y s }3 ± = 0, so that 

(3.29) 9t ± = r lt± x + r 0l± , 

(3.30) y±=p±- 

Lemma 3.1. Let us assume b(, + and b$ + q~ 112 , q 112 . Then the off-diagonal components of the deformation 
matrix are given by 

(3.31) p+ = -a n ri _ + «„?%+, 

(3.32) p_ = -a n ri _ + d n r h+ . 
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Proof. We resolve the A-B compatibility relation ( 12. 52b into monomials of x. Examining the x 7 coefficient 
of the (1, 2) and (2, 1) components yields (13.31ft and (13.321 1 respectively. □ 

Lemma 3.2. Let us assume fa 4- q 1/2 , a„,a n + and A„ + fat, b~^t. Then the spectral and deformation matrices 
satisfy the following residue formulae 



(3.33) 

(3.34) 

and 
(3.35) 

(3.36) 



%+(faqt, qt) 



^{b~ l qt,t) 



W+jbfqt^t) 
Z+(b~ l qt, qt) 



^ + {bl l qt,t) = 0, 



K + (b 6 qt, t) + ^f/'% + (b 6 qt, t) = 0, 
* + (b- l qt,t) + - 6 - jp+gSV,*) - 0. 

-l+(D, 1,1) 



Proof. In this step we compute the residues of the A-B compatibility relation, with respect to x, at the 
zeros and poles of 

(x - qfat) (fax - qt) 



(3.37) 



X(x, t) = 



q(x- fat) (fax - f) 

From the residue of ( 12.53b at the zero x = faqt we deduce 113.331 1, and from the same equation at the zero 
x — b~ x qt we deduce ( 13.34b . From the residue of (12. 56b at the pole x = fat we deduce (13.35b , and from the 
same equation at the pole x = b~ 1 t we deduce (13.36b . □ 

We introduce our first change of variables, u n , v n H> z+, via the relations 
(3.38) 
(3.39) 



V„ = l K + Z + + K - z -\, 



lln = 



[K+Z+ - K-Z-] . 



2(q-l) 

The new variables satisfy an identity corresponding to A3. 141) which reads 

(3.40) k + k-z + z- = — [w 2 - Ay 2 V 2 ] (A n , t). 

A n 

Next we subtract (|3.34b from (|3.33b , in order to eliminate both z_ and Tq-. This yields 



(3.41) 



(l - qbfj a n 



fa(faqtA n - 1) 
fa An 



qb\t 



(%„ - faqt) (b 6 A„ - qtj 



+ qt—n,- = 0. 

fa 



This result motivates the definition of the new variable Z 



(3.42) 



z = 



fafaqt 



-z+ + 



b 5 qtA„ - 1 



(A„ - faqt)(faA„ - qt) 

Definition 3.1. In terms of this new variable Z we have 

1 {A n - tfa)(faA„ - t) [(b 5 t - z)A n - 1] 



(3.43) 
(3.44) 



z+ = 



fafat A„ 
(hA n - l)(b 2 A„ - l)(M« - l)(M» - 1) 



fat 



A n [(fat - Z)A„ - 1] 
Our final rewrite of the dependent variables is 

(3.45) A n (t) ^ g(t), 

(3.46) Z(t)^fat-f(q-% 

We are now in a position to undertake the third stage of our derivation. The first of the evolution 
equations is given in the following result. 
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Proposition 3.2. Let us assume that q + 0, + q 1/2 , t + 0, g + Q,b(,t,b^t and a„ + 0. The variables f ,g 
satisfy the first order q-difference equation 

(g-ht)(g-b-H) 

This evolution equation is identical to the second equation ofEq. (4.15) of Kajiwara et al ]~~] and to the second 
equation ofEq. (3.23) of Kajiwara et al, (6). 

Proof. Subtract l|3.36b from l|3.35b in order to eliminate both z+ and ro,+. This yields the relation 



(3.48) ' hh 



(-l + qbfja n 



t (h - tA„) 

:Z- - 



(b 6 t - A n ) (-f + b 6 A„) " b 6 A„ 



1 

+ qt—n A = 0. 
be, 



Thus we have two different ways of computing the ratio of /%+ to r\-) on the one hand we have from 

CHE) 

\b 5 (qb 5 t- z)-t]a n 

(3.49) r li+ = - L v ' — i—n,-, 

b 5 Za„ 

whereas using <|3.481 > we have 

(3.50) h J5lL r 2± = \ hh f {b\A n - 1) {b z A n - 1) (h A n - 1) (b 4 A n - 1) 

+ (b 5 - tA„) (A„ - tb 6 ) (b 6 A n - t) [A„ (b 5 t - Z) - 1] ] 
- [t z b 6 (hA n - 1) (b 2 A n - 1) (b 3 A n - 1) (M„ - 1) 

- (qb 5 tA n - 1) {A n - tb 6 ) (hA n - t) [A n (b 5 t -Z)-l] ]. 
Equating these two forms gives (13. 47b . □ 
The second evolution equation, to be paired with the first (13.4711 as a coupled system, is given next. 

Proposition 3.3. Let us make the following assumptions: t + 0, £>5 + 1, q~ 112 , q~ l , f + 0, bsf + t, g + 0, g + 
and a„ + 0. In addition let us assume that the condition 

(3 51) . + 1 - qb 5 tg - qb\ + qbjfg 

f - hqt 

holds. The variables f, g satisfy the first order q-difference equation 



(3.52) (f&-l)(fg-l) = qi 



2 (f-b 1 )(f-b 2 )(f-b 3 )(f-h) 
(f-b 5 qt)(f-b-H) 



This evolution equation is the same as the first equation ofEq. (4.15) in Kajiwara et al [5j and the first equation 
ofEq. (3.23) in Kajiwara et al |6j, both subject to typographical corrections. 

Proof. Cross multiplying the relations I l3.33ll3~34ll3.351l3.36b we can eliminate all reference to the defor- 
mation matrix and deduce the identity 

W + (b 6 qt, qt) W_(b 6 t, t) _ W + W?t, qt) W-(b?t, t) 

Z + (b 6 qt,qt) Z + (b 6 t,t) % + {qbft,qt) % + (bft,t)' 

Into this identity we employ the following evaluations for the advanced and retarded values of z± 

f 1 (fg-l){g-hqt){gh-qt) 



(3.54) 

(3.55) z_ = qb 5 t 

(3.56) z 4 
(3.57) 



b 5 b 6 t g 
(gh-l)(gb 2 -l)(gb 3 - 1)^-1) 

g(fg-V 

t (g^-l)(gb 2 -l)(gb 3 -l)(gfe 4 -l) 
'5 g(fg-D 
_ h (fg-l)(g-b 6 t)(gb 6 -t) 
ht g 
n 



We find that this relation f actorises into two non-trivial factors, the first of which is proportional to 

(3 5g) . _ 1 ~ qhtg - qb\ + qbjfg 

f - qfet 

Assuming this is non-zero our evolution equation is then the remaining factor of d3.53t 



2 (f-h)(f-b 2 )(f-b 3 )(f-h) 
(f-qhtXhf-t) 



(3.59) (fg-l)(fg-l) = qb 5 t 

or alternatively J3.52I I. □ 

Lastly we have an auxiliary evolution equation which controls the normalisation of the orthogonal 
polynomial system. 

Proposition 3.4. Let us assume b(, +§,b$ + q~ 112 , f + b$qt, b~ l t and y„ + Ofor n > 0. The leading coefficient 
of the polynomials or second-kind solutions (see i2.13l\2.14$ ) satisfy the first order q-difference equation 

2 f-^t 



(3.60) 



Yn 



^hynj f~b 5 qt' 

Proof. Using the leading order, i.e., the [x] terms, in the expansions d2.38l2.40p with definitions i2.46l2.47t 
we can compute r\ t +. However by using these same expansions to compute /%_ and the equation | |3.49|| , 
which relates these two quantities, we have an alternative expression for Equating these expressions 
then gives d3.60b . □ 

We conclude our discussion by summarising our results for the spectral and deformation matrices in 
terms of the /„, g n variables. Henceforth we will restore the index n on all our variables. The form of the 
spectral matrix is given in the following proposition. 

Proposition 3.5. Assume that \q\ + 1, b^ + 0, b 2 + q~ x , 1, q and a n + 0. The spectral matrix elements ( 12.261 
I2,28ll2.29l > are given by 



(3.61) 

(3.62) 

and 
(3.63) 



b 

= y (q' 1 - bfj a n x (x - g n ) , 
%i- = -r- (l - q _1 bl) a n x (x - g n -i) , 



x(x - g„) 



= -xb 5 b 6 + 



(1 - bi) 



5' I 
6t 



b 5 , / 1 1 1 1 \ 

-T + Mr + r + r + r 

t \b\ b 2 b 3 b±) 



b 6 (b 5 f n -t) 
(l-b 2 )t 



gn(t-f„b 5 ) tb 5 (l + b 2 6 ) 



fn 



(!- fo 5) L 



(1 - gnbl)(l - gnbl){l ~ gnb 3 )(l ~ gnh){gn _ xb\) 
(1 - fngn)gl(x ~ gn) 



b\ 1-bl 



or 



(3.64) 



2B«- 



gn X gn 

b 6 t (fn - h){fn ~ b 2 )(fn - b 3 )(f„ - b A )(l - b\f n x) 



n Jn 



x(x - g„) (1 - bl) 



f$0- ~ fngn)0- - fnX) 



| b 6 t(l - xh)(l - xb 2 )(l - xb 3 ){\ - xb 4 ) b 6 (b 5 f n - t) 



X(l -Xf„)(x-g n ) 



fn 



b b {hf n -t) ( bl(l + b 2 6 )t fr 5 g h( ,b 5 \ /l 1 1 1\ 

TT-, TtTT i Z + — T~ ( ^ ~ b 5Jn) + 75 t + b 5fn-tf„\ — + — + — + — \ 

b 5 {l-bl)t \ h f n fit V&i b 2 b 3 bj 



and 



(3.65) 



f2B«, 



(1 - xf n )(x - b 6 t)(xb 6 - 1) b e (t - hb 5 )(t - b 2 b 5 )(t - b 3 b 5 )(t - b A b 5 ) 



X(X - g„) 
| h{b 5 f n -t) 



X(X - g n ) 



1-b 2 



b 5 (tg„ - h) 



n _ , 2 , _t 2 _ 2 gfl + b\)t _ b 5 t 2 (l - g n h)(l - g n b 2 ){\ - g n b 3 )(l - g„h) 



gn0--fngn)(tgn-fe) 
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or 



f3B„,_ (1 -*/„)(* -fe 6 f)(&6*-f) ifn-h){fn-b 2 )(fn-h){f n -h) 
(3-66) — r = ; ; H 



x(x - g n ) 
| hjbsfu-t) 



(1 - + + — (f - & 5 /„) + — 

D 6 Jn /„ 



/ 1 1 1 1 \ 

t+b5f "- t Hv 1 + v 2 + h + vJ 



Proof. This follows from applying the transformations 113.3813. 39b , (13.4313. 44b and Il3.45l3.46ll successively 
to (13.13b , Il3.11b and Il3.12b . The alternative forms arise from applying partial fraction expansions with 
respect to either of /„ or g„ . a 

The deformation matrix is summarised in the next result. 
Proposition 3.6. Assume that \q\ + \,b(,+ 0, b 2 + q~ x , l,q. The deformation matrix elements are given by 

q (eb 6 + gibib 6 -tg n bi(i + bi)) 



(3.67) %,, + = j^-\x + 



hy n I 1 - b\ 



gnh 

qt 2 b 5 (l - hg n )(l - b 2 g„)(l - hg„)(l - b 4 g n ) 
gnifngn ~ l)(tgn - h) 

q(t - hb 5 )(t - b 2 b 5 )(t - b 3 b 5 ){t - b±b 5 ) 



(tg„ - b 5 )(f n b 5 - t) 



and 



(3.68) 9t„,. 



fn-1 



1-bl 



q [t 2 h + glblb 6 -tg n (l + bl)) 

gnh 



qt 2 b 5 (l - b lgn ){\ - b 2 g n ){l - hg n ){\ - b ign ) 

gnifngn -Wgn-h) 

q(t - hb 5 )(t - b 2 b 5 )(t - b 3 b 5 )(t - b 4 b 5 ) 



(tg n - b 5 )(f n b 5 - t) 



Furthermore 
(3.69) 

(3.70) 





fn 


hyn ' 


a„ 


bbjn 


fn 




fn-1 


hy n -i 


hy n -i 


fn-1 



Proof. Using the leading orders in the expansions (12.3812. 40b , i.e., the [x] terms, with definitions (12.4612. 47b 
we deduce 



(3.71) 
(3.72) 



fn 

hy„ ' 
hy n -i 

fn-1 



Using the leading orders in the expansions (12.39b , i.e., the [x°] terms, with definition l|2.48b and l|2.36b we 
deduce 



(3.73) 
(3.74) 



P+ = a„ 



fn 


hy n 


hy n 


fn _ 


fn-1 


hy n -i 


hyn-i 


fn-1 
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Using the coefficient of the [x 7 ] term in the (1, 1) element of the A-B compatibility relations d2.50b , along 
with the solution of ([3.47b for f n {q~ l t) and l|3.52b for g„(qt) we deduce 



(3.75) 



= qh(y- h) gn + - l) gn + [l + qb 2 5 - h{b~ l + b~ x + b" 1 + fc" 1 )^] j + ^^b\qt, 

whereas if we examine the [x 7 ] term in the (2, 2) element of the A-B compatibility relations in the same 
way then we find 

(3.76) (l-b 2 /^ 
<x- 

= qh L - j\ g n + (l - ^) gn + [b5iK l + K + K l + K'W -qbl-l]j- 

Into both of these expressions we can employ (|3.52b for g„ and make a partial fraction expansion with 
respect to /„. □ 

4. Reconciliation with the Lax pairs of Sakai and Yamada 

4.1. Sakai Lax Pair. In (TBI Sakai constructed a Lax pair for the Eg q-Painleve equations using a 

degeneration of a two-variable case of the Gamier system based upon the Lax pairs for the D 5 q- 
Painleve system (141 . Subsequently Murata [10] gave more details for this Lax pair. We intend to 
establish a correspondence between our Lax pair and that of Sakai. We will carry this out in a sequence 
of simple steps rather than as a single step as this will reveal how similar they are. 

Our first step is to give a variation on the parameterisation of the spectral and deformation matrices 
to that given in Section [3] In this alternative formulation, we seek a spectral matrix A(x;t) (actually 
identical to the Cayley transform A* n ) with the specifications 

(4.1) A(x; t)=A + Aix + A 2 x z + A 3 x 3 , 

and 

(i) The determinant is 

b(, {b\x - 1) (b 2 x — 1) (b 3 x - 1) (b^x - 1) (x - b^t) (b^x - t) , 

(ii) A3 is diagonal with entries K\ = -b^b^ and k 2 = —b(,/b$, 

(iii) A = b 6 tl, 

(iv) The root of the (1,2) entry of A(x; t) with respect to x is A, 

(v) A(A;f) is lower triangular with diagonal entries -b^b(,Az+ and — ^ — where b 2 Z-Z + A 2 = 

"5 

detA(A; t). 
Any such matrix is in the general form 



(4.2) A(x; t) = tb 6 I - 



(b 5 b 6 x \z\ + (x - a) (x - A)] b 6 wx(x — A) 

b 5 b(,x(xy + 5) b^x [z 2 + (x- f$)(x- A)] 

w b~ 5 



where the properties specify the variables 

(4 .3) (1 _^ = l + l + l + l + (l + , 6 ) t 

t 1 1 \ 1 b\z\ z , 
- (h + b 2 + b 3 + b 4 )b 5 j -b 5 y— + b 6 j - + -j- + j-2A, 

(4.4) (i-b^=-(i i+ i 2+ i + i-)bi-{i 6+ b 6 y 5 t 

t 1 1 \ 1 b\z\ z, 
+ (h + b 2 + b 3 + h)b 5 j + b 5 \j- + b 6 j - - -A- - — + 2b 2 5 A, 
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(4.5) y = - (b 3 b A + b 2 b 3 + b 2 h + fab 2 + hb 3 + b x b±) - + — + — + — j Q- + fc 6 j ( 

+ ajS + zi + z 2 + 2(a + /3)A + A 2 , 



(b 3 b 4 + &2&3 + fafa + M2 + fab 3 + b t h) (— + £> 6 ) ~{t- + fa) t 
1 1 1 1 \ 

— + — + — + — U 2 - Zi(j8 + A) - z 2 (a + A) + (-2ajS + y)A - (a + |3)A 2 
Oi b 2 b 3 bi) 



r 2 ,i r 2/2 



(4.6) 5 = b\ + b 2 + b 3 + bi + 

+ 

The z\ and z 2 are related to z± by 

(4.7) zi = z+ + — ^— , and z 2 = z_ + — . 

05A A 

In addition 

\-qb\ 

(4.8) a; = — a n . 

<\ 

We seek a deformation matrix B(x; t) of the form 

(4.9) B(x;t) = - — £ ——(xl + Bo), where B = 

(x - beqt)(x - b 6 L qt) 

This leads to the compatibility relation 

(4.10) B(qx; t)A(x; t) = A(x; qt)B{x; t). 

This relation is just a rewriting of (12.50b whereby all the factors of x are placed into the denominator of 
S by the above definition. 

Lemma 4.1. The overdetermined system (14.101 , with d4.1|l and (14.9b zs satisfied if the coupled q-Painleve 
equations i3A7i and 1 13.5211 are satisfied. 

Proof. Examining the coefficient of x 6 in the numerator of the (1,2) entry of (14. 10b we find 

(4.11) r u = -_i_(t2>- ro ). 

Now we seek two alternative expressions for r\ j2 — one involving quantities at the advanced time qt 
and another involving those at the unshifted time t. The first of these is found from solving for the (1,2) 
entry of the residue of (|4.10b at x = b^qt simultaneously with the (1,2) entry of the residue of (|4.10b at 
x = b~ l qt. This yields 

-qtfv(qtb 6 - A)(qt - b 6 X) 

(4.12) r 1/2 - 



b 5 \qt(l - fr 5 Mi) + (#6 - A) [fa + qtfa(qt - M)]} 

The second expression for r\ i2 is found from solving for the (1,2) entry of the residue of (14.10b at x = b^t 
simultaneously with the (1,2) entry of the residue of l|4.10b at x = b~ l t. This gives 

r = -gtw(b 6 t - A)(t - M) 

j ri ' 2 tb 5 -tb 6 z 2 + b 5 b 6 (b 6 t-A) + t(b 6 t-A)(t-b 6 A)' 

Combining d4.11|l and (|4.12b or ll4.11b and (14.13b , and employing the change of variables d3-43b and (13.44b 
with (|3.45b and (|3.46b , we can solve for w in two ways. Assuming w is non-zero it cancels out, leaving 
an expression for Z in terms of Z and A. This is equivalent to the first ^-Painleve equation l !3.47b . 
To find the second equation we solve d4.10b for A(x; t) 

(4.14) A(x; t) = B(qx; t)~ 1 A(x; qt)B(x; t), 

and use this to find the zero of A(x; t)i 2 , i.e., g(t). In addition to r\ r2 we now require r 2/2 (even though the 
denominator of A(x; t)i 2 depends on r\ r \, r\ i2 , r 2i \, r 2j2 identities resulting from the compatibility conditions 
imply that this will trivialise - see the subsequent observation). The entry r 2j2 has already been found, 
along with ri i2 , from the arguments given earlier and this is 

^ 5 ^ 1 + h D n, t JV( fc3 , b 6 + qifb 5 (l + bl)-b 5 b 6 qt(a + A) 

(4.15) r 2 ,2 = A (b 6 qt - A)(qt - b 6 A)- 



bj(A - faqt) - b 6 qt + fab\qtz x - b 5 b 6 qt{b 6 qt - A)(qt - b 6 A) 
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The numerator of A(x; t)u appears to be a polynomial of degree 6 in x, however it has trivial zeros 
matching those of the denominator 

q 2 b$ 

— — (x - b 6 t)(x - b 6 qt)(b 6 x - t)(b 6 x - qt)w, 
b l 

so that their ratio is in fact polynomial of degree 2. Into A(x; t)u we first substitute for r\£ using l!4.11b , 
then for ^2 using (14.151 , and thirdly for &, /?, y, 5 using J4-3U4.4ll4.5U4.6b at up-shifted times, respectively. 
Into the resulting expression we employ ii.7\ for z\,Z2 along with ||3.43I3.44| | at the up-shifted time to 
bring the whole expression in terms of A and Z. The relevant zero of the ensuing expression (the other 
zero is x — 0) then gives A = g in terms of A = g and /, or equivalently by (13.52b - □ 

Now we recount the formulation given by Sakai [15] and Murata [10J. Their Lax pairs are 

Y{qx; t) = A(x, t)Y{x; t), 
Y(x; qt) = B(x, t)Y(x; t), 

satisfying the compatibility condition A(x, qt)B(x, t) = B(qx, t)A(x, t) . The spectral matrix is parameterised 
in the following way 

subject to the key properties 

(i) The determinant of A(x, f) is 

Ki%2(x - a\){x - az){x - a^)(x - a.i)(x - a^,t){x — a^t), 

(ii) A3 is diagonal with entries Ki and K2 = qK\, 

(iii) Aq has eigenvalues 9\t and 02t, 

(iv) The single root of the (1,2) entry of A(x, t) in x is A, 

(v) A(A, t) is lower triangular with diagonal entries K1/J1 and K2f/2- 
Given these requirements, the entries of A(x, t) are specified by 

L(x, t) = x - A, 

Z(x, t) = \ii + (x- A) [62 + x 2 + x(y + A)] , 
W(x, t) = [i\ + (x - A) \6i + x 2 + x (-y -e 1 + A)] , 
X(x, t) = [WZ — (x — a\) (x - ai) (x — 03) (x - 84) (x — a^t) (x — a^t)] L -1 , 



where 



(ki - K2)6i = A 1 [KijLii + K2/i2 — B\t — 6ii\ - K2 \y {y + £i) + 2A 2 - Ae\ + e 2 ] , 

(k\ - Ki)b2 = -A -1 [xij.ii + K2H2 - 0\t - d 2 t] + ki \y (y + e\) + 2A 2 - Ae\ + £2] / 
f/ifi 2 = (A - fli) (A - a 2 ) (A - fl 3 ) (A - fl 4 ) (A - a 5 t) (A - a 6 t) , 

Here e, is the / h elementary symmetric function of the indeterminates [ai, 0.2, as, a^, ast, a^t] . Despite the 
expression for X(x, t), it is a quadratic polynomial in x. In Murata's notation we have p. = \i = fj-2, 
5 = 61 and 6 = 62- The deformation matrix B(x, t) is a rational function in x of the form 

B(x t) = x(xl + Bo) 

(x - a 5 qt)(x - a 6 qt) 

Next we consider the first transformation of the Sakai linear problem with the following definition: 

y(x,t) = s(x r ty 1 Y(x,t), 

and 

' = { I ° 

\S\ + S2X x 

The transformed spectral linear problem is 

y{qx, t) = a(x, t)y{x, t), 
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with a transformed spectral matrix 

A{x, t ) = S(qx, t)~ 1 A(x, t)S(x, t). 

We fix the parameters of the transformation by the requirement that the coefficient of x~ l in the (2,1) 
entry of S is zero (only the (2,1) entry is non-zero) and also that the coefficient of x° in the (2,1) entry of 



S is zero. Thus we find 



1 

si = 2cjKlZV/ { W 02 ~ e ^> t + K i 0*1 " W'2 + (^2 - 6i)A)] , 



(20 \ 
— —t - c\0\\\2 - 02^i j A" 2 - qx\e^t~ x A~ x -e\0 2 + {q0\ - 0i)y + {q0\ + di)A. 



120102, 
Ki 

The new spectral matrix can be parameterised by the polynomial 

JA(X, t) - + AlX + A 2 X 2 + J?h,X 3 , 

and possesses the following properties 
(i) The determinant of J{(x, t) is 

k 2 (x - fli)(x - ci2){x - a$)(x - Oi)(x - a^t){x — a^t), 

(ii) J4. 3 = Kit, 

(iii) J%q is diagonal with entries 0\t and q~ 1 02t, 

(iv) The roots of the (1,2) entry of A in x are and A, 

(v) A(A, £) is lower triangular with diagonal entries K\[i\ and K2/U2. 
Any such matrix admits the general form 

. _ l6\t + K\X [(x - A) (x - a) + Vl] ^Kiwx(x-A) 
^' ' ^ 7CiZf _1 x(xc + (£) <? _1 02f + Kix[(x - A) (x - 6) + v 2 ]j' 

where the properties given above fix the introduced parameters as 

{q0\ - 02)a = [02Vi + q0\V2 + ^K^r 1 ] A -1 + q0\e\ - 2q0\A, 

{q0\ - 2 )6 = - \02V\ + q0\v 2 + ^K^r 1 ] A" 1 - 2 e\ + 20 2 A, 

qc = ab + 2(a + £)A + A 2 - £2 + + V2, 

^ = -(a + 6) A - 2a6A - aV2 - bv\ + (qc — v\- V2M + £3 H f. 



The variables, vi and V2 are defined by 

Kl/ij - 01 1 (JK1/J2 - 02^ 

vi = , and v 2 = - . 

KiA ^KiA 

The transformed deformation matrix H is computed using 

®(x, t) = S(x, qt)~ x B(x, t)S(x, t), 

and has the form 

X(x3n + 1) 

® = 



(x - a 5 qt)(x - a b qt) 
We define a new variable v using 

(A - «i)(A - a 2 )(A - a 3 )(A - 04) 

f2 = j — ; , 

A - v 

and by implication 

III = (A - a 5 0(A - fl 6 t)(A - v). 

Using identical techniques to those employed in the proof of Lemma |4~T1 we can show that the compat- 
ibility relation leads to the evolution equations 

(A - «i)(A - a 2 )(A - « 3 )(A - 04) 



(4.16) (A - v)(A - v) 

<4 - 17) ('-jIH)- 



(A - a 5 f)(A - M) 
Xj\ A) q (a 5 a b tv + 0ilqKi)(a 5 a b tv + 2 lqKi)' 
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To make the full correspondence between our system and this one we must consider a further 
transformation, given by the linear solution 

The prefactors are elliptic functions defined in terms of the ^-factorial by 

& q (z) = [q,-qz,-z ;q) oo , e,, t (z) = _^ , 

with properties 

9 q (qz) = qz$ q (z), e q>i (qz) = te q , t (z), e Mf (z) = ze ¥ (z). 
This is the solution satisfying the linear equations 

V(q~ 1 x,t) = H(x,tmx,t), 
y{x l q- 1 t) = to{x,tMx,t). 
The transformed spectral matrix s 2t is given by 

(4.18) 3t(x, = fx^x -1 , r 1 ) = 3I 3 + % 2 x + «jx 2 + % x 3 

which swaps the roles of the leading matrices around x - and x = oo. This spectral matrix has the 
properties 

(i) The determinant of ?l(x, t) is 

7Cj(l - fliX)(l - fl2X)(l - fl3X)(l - «4X)(f - fl5X)(f - flgX), 

(ii) 5l 3 = Klfl, 

(iii) 3Io is diagonal with entries 0j and 02/ 

(iv) the roots of the (1,2) entry of *2l(x, t) in x are and A -1 , 

(v) 9t(A, f) is lower triangular with diagonal entries Ki/./if A" 3 and Ki/^fA -3 . 

The transformed deformation matrix has the form 

Since the compatibility relation between 1)4.181 and l|4.19t is rationally equivalent to that for f, the 
evolution equations are the same. 

It is clear that s 2) and Y n satisfy equivalent linear problems and that the following correspondences 
hold: 

q i— > q~ l , 
t h-> r\ 

v(t) k» /(r 1 ), 

l 

Ki i-> 06/ «i ht (i = 1/2,3,4), a 5 i-> — , 
0i i-> -e 5 »6/ q~ x 6z^>—r-- 

4.2. Reconciliation with the Lax pair of Yamada fTTI . In his derivation of a Lax pair for the Eg q- 

Painleve system Yamada employed the degeneration limits of Eg 1 ' ^-Painleve — > (j-Painleve — > Eg 1 ' 
q-Painleve. In doing so he retained eight parameters b\, . . ., bg constrained by qb\bjb^b\ = b^b^ojbg, and 
his E^ 1 ' ^-Painleve equation was given by the mapping of the variables 

t^q'H, f,g^f,g, 
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subject to the coupled first-order system (see his (36)) 
(4.20) 



if 8 ~ Wg ~ i) (kg ~ Vihg ~ Vihg - i)(hg - 1) 

-=q— 



(4.21) 



// ' hhihg - t)(b sg - t) 

ifg - Wg - 1) (h - f)(b 2 - f)(h - f)(h - f) 



gg 



(j-b 5 t)(f-b 6 t) 



The Lax pairs constructed by the degeneration limits were given as a coupled second-order ^-difference 
equation in a scalar variable Y(z, t) (see his (37)) 



(4.22) 



(hq - z)(b 2 q - z)(b 3 q - z)(hq - z)i 
q(qf - z)z 4 



2 r 



Yiq-h) 



g- 



fi(gz-q) 



Y(z) 



q(hg - l)(b 2 g - l)(hg - l)Q>tg - 1) b 5 b 6 (b 7 g - t)(b s g - t) 



g(fg - l)z 2 (£z - q) 



and a second-order, mixed ^-difference equation, 



fgz 3 

(b 5 t - z)(b 6 t - z) 
t 2 z 2 (f - z) 



Y(qz) 



Y(z) 

t z (gz-l) 

g Z 



Y(z) 



= 0, 



(4.23) 



gz 
t 2 



Y(z) + (q- gz)Y{q- l z) - q~ 2 gz{qf - z)Y{q~ l z) = 0. 



In order to bring (|4.20l l and (|4.2H into correspondence with our form of the E 6 ^-Painleve system 
(see dl.H and dl.2|l ) we will employ the following transformation of Yamada's variables 



t ^ r 1 , 



zhz , 



f,g^g~\r\ m = Y( Z - i ) / 

and the specialisations of the parameters 

b 5 i-> b' 1 , b 6 i-> b 6 , b 7 m> qb 5/ b 8 i-> b~ l , 

so that b$b(, = 1 and b 7 b% = q. Under these transformations we deduce that (14.20b becomes l ll.2|l and 
d4.2H becomes {1.1) . Furthermore the pure second-order divided-difference equation l|4.22|l becomes 

nLd - bjz) . 

(4.24) —L- — Y(z) 



t 2 z(z-g) 

U%i(l -bp) zU%x(bj-f) z(f-b 5 qt)(b 5 f-t) (z - b 6 qt)(b 6 z - qt)(q -fz) 



z(z-g)(l-/z) /(l-/g)(l-/z) 



hqt 2 f 



b 6 qt 2 z(z - qg) 



Y{q-h) 



(z - b 6 qt)(b 6 z - at) _ , 

+ i / v y(q~ 2 z) = o, 

b 6 z{z - qg) 
Y(z;qt) = 0. 



and the mixed divided-difference equation (14.23ft becomes 

qt 2 . . (1 - fz) . (z - <?) 

(4.25) ^(f V) - q K -j^Y{z;t) - ^ 

Having put Yamada's Lax pairs into a suitable form we now seek to make a correspondence with our 
own theory and results. A single mixed divided-difference equation can be constructed from the matrix 
Lax pairs ( (12.171 and d2.35l l). For generic semi-classical systems on a ^-lattice grid we can deduce either 



11 1 

~W + AyV<p7 P "^ ; ^ + (W + AyV)(R + AuS) 
or an alternative, 

Z+(x) 



p„(x; f) + 



1 



R + AuS 31 



p„(qx;t) = 0, 



(4.26) 



(W-AyV)(x) 
+ 



p„(qx; qt) 



1 



(W - AyV)(x)(R + AuS)(qx) 



[X + (x)K + (qx) + y + (qx)W-(x)]p n (qx;t) 

%K(qx) 



(R + AuS)(qx) 



p„{x;t) = 0, 
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which we will work with. Using the spectral and deformation data d3.6l3.25| l and the explicit evaluations 
of the deformation matrix J3.67II and spectral matrix Q.65II , we compute the coefficients of the above 
equation 

l-qbl 

(4.27) -(R + AuS)(qx)Z + (x) = a n — -^-x(x - b 6 t)(x - g), 

05 

(4.28) -(W - AyV)(x)¥ + (*0 = ^ " ^ ^ " ^ " M( * " M) 



(4.29) Z + (x)ft + (qx) + ^ + (qx)3B-(x) = 



fn h faqt-f 

a„y„ b 6 (l ~ qb 2 5 ) (x - b 6 t)(b 6 x - f)(l - fx) 



fn h b 5 qt - f 



Now we set p n = FU where F is a gauge factor and U is the new independent variable, into (14.26b and 
make a direct comparison with (|4.25t . Comparing the coefficients of U(x; t) and U(qx; t) in this later 
equation we deduce that 

r4 , m Hqx.t) _ 1 (l-hxxi-hx) 

( 1 F(x,t) f 2 1-hxt-^ ' 

A solution is given by 

(b ( ,xr 1 ;q) 
' (b 2 x, b 3 x; q) D 



(4.31) F(x, t) = e qit -2(x) , u z u C(x, t), 



where C is a ^-constant function, C(qx, t) = C(x, t).Now comparing the coefficients of U(qx; qt) and U(qx; t) 
in the previous equation we find that 

(4 32) F(qx,qt) _ y„ b 6 (b 6 x - t) 

F(qx, t) f„ qg(b 5 qt - f)x 2 ' 

Substituting our solution into this equation we find a complete cancellation of all the x dependent factors 
resulting in a pure ^-difference equation in t 

(4.33) f " C - hqt 



7nC g(f-b 5 qt)' 

Thus we just need a solution C(f ) independent of x, however we only require the existence of a non-zero, 
bounded solution rather than knowledge of a specific solution. In conclusion we find that our new 
mixed, divided-difference equation is now 

(4.34) t 2 U(x; t) - (1 - fx)U(qx; t) - ?—^-U(qx; qt) = 0, 

qgx 

which is clearly proportional to J4.25II with the identification U(x; f) = Y(q~ 1 x; t). 

A second-order (^-difference equation in the spectral variable x for one of the components, say p„, was 
given in J2.30|l , and for ^-linear grids can be simplified as 

W + AyV r2B. 3DB_ , 1 W - AyV , 

(4.35) * (x)p n (qx) - -±(jc) + -=-{jf x x) - '-^ y 



p„(x) + — = — (q lx )Vn{q x x) = 0. 



From the explicit solution of the gauge factor we note 

F(qx,t) (1 - b 2 x)(l - b 3 x) Htf x x,t) t{t-b(,q- l x) 



(4.36) 



F(x,t) t(t-b 6 x) ' F(x,t) (l-b 2 q- l x)(l-hq- 1 x)' 

Substituting the change of variables into j4.35|l we compute that 



W + AyV F(qx, t) 1 b 6 Yl%iQ- ~ tyx) 



(4.37) - (x) 



F(x,t) (q-l)uia„ t x(x-g) 



W-AyV _j F{q- l x,t) _ 1 t(x - b 6 qt)(b 6 x - qt) 
( ' % + W X) F(x,t) (q-l) Ul a„ x(x-qg) 
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In addition, using the explicit representations of the diagonal elements of A*, i.e., 3B± (see (13.63b , 113, 65b ) 
we compute that 



(439) ~ hi 



3B 4 



3B_ 



x(x - g) x(x - g) 



U.%i(l-bjX) xU%i(f-bj) x(f-b 5 qt)(b 5 f-t) (x-b 6 qt)(b 6 x-qt)(q-fx) 



x(x-g)(l-fx) f(l-fg)(l-f X ) 
In summary we find 



b 6 qt 2 x{x - qg) 



(4.40) 



U(qx) 



U%iO- - bjx) 
t 2 x(x - g) 

U%i0--bjx) xll%i(f-bj) x(f-b 5 qt)(b 5 f-t) (x-b 6 qt)(b 6 x-qt)(q-fx) 



x(x-g)(l-fx) /(l -/*)(!-/*) 



hqfif 



U(x) 



U{q~ l x) = 0. 



b 6 qt 2 x(x - qg) 

(x - b 6 qt)(fox - qt) 
b 6 x(x - qg) 

Thus we can see that (14.40b agrees exactly with (|4.24b and the identification noted above. 
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